Scaling laws of creep rupture of fiber bundles 



(N 
O 

o 

(N 
D 

cn 



Ferenc Kun^*, Raul Cruz Hidalgo'^, Hans. J. Herrmann'^, and Karoly F. PaP 
^Department of Theoretical Physics, University of Debrecen, 
P.O.Box: 5, H-4010 Debrecen, Hungary 
Institute for Computational Physics, University of Stuttgart, Pfaffenwaldring 27, 70569 Stuttgart, Germany 
'^Institute of Nuclear Research (ATOMKI), P. O. Box 51, H-4OOI Debrecen, Hungary 

(February 1, 2008) 

We study the creep rupture of fiber composites in tiie framework of fiber bundle models. Two 
novel fiber bundle models are introduced based on different microscopic mechanisms responsible for 
the macroscopic creep behavior. Analytical and numerical calculations show that above a critical 
load the deformation of the creeping system monotonically increases in time resulting in global 
failure at a finite time t j , while below the critical load the system suffers only partial failure and the 
deformation tends to a constant value giving rise to an infinite lifetime. It is found that approaching 
the critical load from below and above the creeping system is characterized by universal power laws 
when the fibers have long range interaction. The lifetime of the composite above the critical point 
has a universal dependence on the system size. 

PACS number(s): 46.35.+Z, 46.50. +a, 62.20.Mk 
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I. INTRODUCTION 



Under high steady stresses fiber composites may un- 
dergo time dependent deformation resulting in failure 
called creep rupture which limits their lifetime, and 
hence, has a high impact on the applicability of these 
materials in construction elements. Both natural fiber 
composites like wood [1-4] and various types of fiber rein- 
forced composites [5-10] show creep rupture phenomena, 
which have attracted continuous theoretical and experi- 
mental interest over the past years. The underlying mi- 
croscopic failure mechanism of creep rupture is very com- 
plex depending on several characteristics of the specific 
types of materials, and is far from being well understood. 
Theoretical studies encounter various challenges: on the 
one hand, applications of fiber composites require the de- 
velopment of analytical and numerical models which are 
able to predict the damage histories of loaded compos- 
ites in terms of the specific parameters of constituents. 
On the other hand, creep rupture, similarly to other rup- 
ture phenomena, presents a very interesting problem for 
statistical physics, it is still an open problem to embed 
creep rupture into the general framework of statistical 
physics and to understand the analogy between rupture 
phenomena and phase transitions. 

Creep failure tests are usually performed under uniax- 
ial tensile loading when the specimen is subjected either 
to a constant load tio or to an increasing load (ramp- 
loading) and the time evolution of the damage process 
is monitored by recording the strain e of the specimen 
and the acoustic signals emitted by microscopic failure 
events [1-10]. In the present paper we study the creep 



rupture of fiber composites in the framework of fiber bun- 
dle models [11-24] focusing on general characteristics of 
creep rupture. Our main goal is to reveal the univer- 
sal aspects of a creeping system which do not depend on 
specific material properties. We will derive scaling laws 
which emerge when approaching the failure stress both 
from below and above, furthermore, the finite size scaling 
of the time to failure is studied. Beyond the theoretical 
understanding, exploring universal features of creep helps 
to evaluate experimental data and can serve as a guide 
to extract the relevant information invariant under the 
variation of specific material properties. 



We consider only the case of global load sharing 
(GLS) for the redistribution of load following fiber failure 
[13,18-23]. For several types of materials GLS provides 
an adequate approach, and it has the advantage that 
many of the important GLS results can be obtained in 
closed analytic forms. During creep at the critical point 
and above it, stress localization occurs which can be cap- 
tured by localized load sharing in fiber bundle models 
[13,17,21] but this is beyond the scope of our investiga- 
tions. The present study is focused on analytic modeling 
in the framework of GLS. We also develop efficient simu- 
lation techniques for testing the analytic results, further- 
more, to study finite size system. Since similar studies 
with local load sharing require large scale computer sim- 
ulations and a large amount of data processing it will be 
presented in a forthcoming publication. 
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II. MODELS 

In order to work out a theoretical description of creep 
failure we introduce two novel fiber bundle models im- 
proving the classical fiber bundle models [11,12] which 
have proven very successful in the study of fracture of 
disordered materials [13-20,23]. The two models differ 
in the microscopic mechanisms assumed to be responsi- 
ble for the macroscopic creep behavior: (i) in the first 
approach the fibers themselves are viscoelastic and they 
break when their deformation exceeds a stochastically 
distributed threshold value (strain controlled breaking), 
(m) in the second model the fibers are linearly elastic un- 
til they break stochastically in a stress controlled way, 
however, after breaking their relaxation is not instanta- 
neous but the sliding of the broken fiber with respect to 
the matrix material, and the yielding or creeping of the 
matrix introduces an intrinsic time scale for the relax- 
ation. 

The mathematical formalism of the first model (vis- 
coelastic bundle) is more simple, hence, the steps of anal- 
ysis will be presented in details for this case. For the sec- 
ond model the corresponding results will be summarized. 



A. Viscoelastic fiber bundle 



1. Analytic model 




FIG. 1. The viscoelastic fiber bundle: intact fibers are 
modeled by Kelvin- Voigt elements. After fiber breaking the 
corresponding element is removed from the model. 



Our model consists of A'' parallel fibers having vis- 
coelastic constitutive behavior. For simplicity, the pure 
viscoelastic behavior of fibers is modeled by a Kelvin- 
Voigt element which consists of a spring and a dashpot 
in parallel (see Fig. 1) and results in the constitutive 
equation 



where /? denotes the damping coefficient, and E the 
Young modulus of fibers, respectively. Eq. (1) provides 
the time dependent deformation s{t) of a fiber at a fixed 
external load (Tq 



1 - e-^'/^ 



-Et/p 



(2) 



where £o denotes the initial strain at t = 0. It can be 
seen that e{t) converges to (lo/E for t oo, which 
implies that the asymptotic strain fulfills Hook's law. 
To incorporate breaking in the model we introduce a 
strain controlled failure criterion for fibers: a fiber fails 
during the time evolution of the system if its strain 
exceeds a damage threshold e^, which is an indepen- 
dent identically distributed random variable of fibers 
with probability density p{ed) and cumulative distribu- 

tion P{£d) — I p{x)dx. Due to the validity of Hook's 
Jo 

law for the asymptotic strain values, the formulation of 
the failure criterion in terms of strain instead of stress im- 
plies that under a certain steady load the same amount 
of damage occurs as in the case of stress controlled fail- 
ure, however, the breaking of fibers is not instantaneous 
but distributed over time. When a fiber fails its load 
has to be redistributed to the intact fibers. Assuming 
global load sharing, the time evolution of the system un- 
der a steady external load CTo is finally described by the 
equation 



1-P{e) 



= (3i + Ee, 



(3) 



where the viscoelastic behavior of fibers is coupled to the 
failure of fibers in a global load sharing framework [26]. 

For the behavior of the solutions of Eq. (3) two distinct 
regimes can be distinguished depending on the value of 
the external load ao'- When do is below a critical value 
ac Eq. (3) has a stationary solution Ss, which can be 
obtained by setting £ = in Eq. (3) 



ao = Ees[l- P{e, 



(4) 



(To = /?£ + Ee, 



(1) 



It means that until this equation can be solved for Ss 
at a given external load ctq, the solution e{t) of Eq. (3) 
converges to Eg when t ^ oo, and no macroscopic failure 
occurs. However, when (Jq exceeds the critical value ac 
no stationary solution exists, furthermore, i remains al- 
ways positive, which implies that for a > ac the strain of 
the system e{t) monotonically increases until the system 
fails globally at a time t / [26] . 

In the regime cTq < ac Eq. (4) also provides the asymp- 
totic constitutive behavior of the fiber bundle which can 
be measured by controlling the external load ctq and let- 
ting the system relax to £«. It follows from the above ar- 
gument that the critical value of the load ac is the static 
fracture strength of the bundle which can be determined 
from Eq. (4) as ac = Esc[i- — P{sc)], where Sc is the so- 
lution of the equation dao/dssl^ = [18]. Since ao{£s) 
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has a maximum of the value ctc at Ec, in the vicinity of 
£c it can be approximated as 



- A{ec - Esf 



(5) 



where the multiplication factor A depends on the proba- 
bility distribution P. A complete description of the sys- 
tem can be obtained by solving the differential equation 
Eq. (3). After separation of variables the integral arises 



t = (3 de 



1 - P{e) 



' ao-Ee[l- P{e) 



(6) 



where the integration constant C is determined by the 
initial condition e{t = 0) = 0. 

The creep rupture of the viscoelastic bundle can be in- 
terpreted so that for CTq < <Tc the system suffers only a 
partial failure which implies an infinite lifetime tf = oo 
and the emergence of a macroscopic stationary state, 
while above the critical load > cr^ global failure occurs 
at a finite time tf, which can be determined by evaluating 
the integral Eq. (6) over the whole domain of definition 
of P(e). 

Below the critical point CTq < <Tc the bundle relaxes 

to the stationary deformation Sg through a decreasing 
breaking activity. To find the characteristic time scale of 
this relaxation process the behavior of s{t) has to be ana- 
lyzed in the vicinity of Eg- It is useful to introduce a new 
variable 6 as 5{t) = — £{t). The governing difi^erential 
equation of 6 can be obtained from Eq. (3) by expanding 
it around 



dt 



E 

'J 



1 - 



£sP{£s) 
L - Pies) 



(7) 



The solution of Eq. (7) has the form S ~ exp[— i/r], 
where r is the characteristic time scale of the relaxation 



process 



E 



1 - P{es) 



It is a very important question how the relaxation time r 
changes when the external driving approaches the critical 
point (Tc from below. Based on Eq. (5) it can be simply 
shown that 



r ^ (C7c — (To) 



-1/2 



for CTo < (Tc 



(9) 



which means that approaching the critical point from be- 
low the relaxation time of the system diverges according 
to a universal power law with an exponent —1/2 inde- 
pendent on the form of disorder distribution. Note that 
a similar power law divergence of the number of succes- 
sive relaxation steps was found in Refs. [27,28] for a dry 
fiber bundle subjected to a constant load. 

Above the critical point the behavior of the lifetime 
of the system can be analyzed analogously when ctq goes 



to (Tc from above. When ctq is in the vicinity of (Tc, i.e. 
(To = (Tc + A(To, where Actq << ac, it can be expected 
that the curve of e{t) falls very close to Sc for a very long 
time and the breaking of the system occurs suddenly. 
Hence, the total time to failure, i.e. the integral in Eq. 
(6), is dominated by the region close to £c when Actq is 
small. Making use of the power series expansion Eq. (5) 
the integral in Eq. (6) can be rewritten as 



tfK. (3 de 



I' 



P(£) 



(10) 



•A<7o-A(£c-£)2' 

which has to be evaluated over a small e interval in the 
vicinity of £c. After performing the integration it follows 

*/ ~ ((To - (Tc)"^/^ for (To > (Tc. (11) 

Thus, tf has a power law divergence at Uc with a uni- 
versal exponent — i independent of the specific form of 
the disorder distribution P(£), similarly to r below the 
critical point. 



2. Finite size effect 

In the above analytic treatment the size of the sys- 
tem, i.e. the number of fibers in the bimdlc, is infinite. 
However, it can be expected that the lifetime of a finite 
bundle has a non-trivial size scaling even in the case of 
global load sharing. During the creep rupture process the 
fibers of a finite bundle break one-by-one in the increasing 
order of their breaking thresholds. Let ej, j = 1, . . . ,N 
denote the breaking thresholds assigned to the fibers in a 
realization of the bundle. Since fibers break one-by-one, 
the actual load on the intact fibers after the failure of i 
fibers is ai = aoN/{N — i) where i = 0, . . . , N — 1, and 
the time At{ei, £i+i) between the breaking of the ith and 
i + 1th fibers (in the ordered series) reads as 



At{si,Si+i) = - 



E 



(12) 



The lifetime tf of a sample of N fibers takes the form 



N-l 

tf{N)=J2Mei,Si+i). 

i=0 



(13) 



In order to determine how tf depends on N, Eq. (13) 
has to be averaged over many realizations of the disorder 
distribution, which can be performed analytically. The 
details of the analytic calculations are summarized in the 
Appendix. Finally, the average lifetime {tf{N)) of a bun- 
dle of N fibers can be cast in the form 

dsEeP{e) [1 - P{e)] 



{tf{N)) « f/(oo) + 



/3(To f deEei 



(14) 



Ee{l-P{e))f 

Eq. (14) shows that for finite bundles the average lifetime 
{tf{N)) converges to the lifetime of the infinite bundle 
tf (oo) as ~ 1/N with increasing number of fibers N. It 
is interesting to note that in the case of global load shar- 
ing the average strength of the bundle ac does not have 
a size dependence. 
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3. Simulation technique 



Most of the analytic results of the previous sections, 
except for the finite size scaling of tf, were obtained for 
infinite bundles. Computer simulations of the creep rup- 
ture of finite bundles are needed to justify the validity of 
analytic predictions for finite systems, and to be able to 
model the rupture process of realistic finite systems. 




obtained by simulations is presented in Fig. 2 for sev- 
eral different values of ctq below and above CTg. The two 
regimes depending on the value of the external load (To 
can be clearly distinguished. 




FIG. 3. The lifetime of a bundle of lO'^ fibers as a function 

of (To — (Tc above the critical point for three different disor- 
der distributions, i. e. uniform and WeibuU distributions with 
p = 2,7 have been considered. The straight line of slope —0.5 
is drawn to guide the eye. 



FIG. 2. e{t) for several values of Oo below and above Gc 
for a bundle of lO'^ fibers. The critical strain and the time 
to failure t/ for one example axe indicated, to denotes the 
characteristic time to = P/E of the system. 



In the framework of GLS an efficient simulation tech- 
nique can be worked out for the failure process. Based 
on the arguments of the previous subsection, the GLS 
simulation of the creep process of a bundle of N fibers 
proceeds as follows: {i) random breaking thresholds 
Ei, i = 1, . . . , N are drawn from a probability distribu- 
tion p, then the thresholds are put into increasing order. 
(m) The time At(ej, ej+i) between the breaking of the 
ith and i + 1th fibers is calculated according to Eq. (12). 
{Hi) The time elapsed till the breaking of the ith fiber 
is obtained as t{ei) = X]}=o from which the 

deformation as a function of time e{t) can be determined 
by inversion. The lifetime of the finite bundle can be 
obtained making use of Eq. (13). 

For simulations we considered a uniform distribution 
of failure thresholds between and a maximum value 
with the probability density function p{e) = l/sm, and 
distribution function P(e) = e/sm- In this case the sta- 
tionary solution, the critical load and the corresponding 
critical strain can be obtained as CTq = Ee[l — s/sm], 
(Jc = Esm/^i £c = £m/2, respectively, = 1 was set in 
all the simulations. The deformation- time diagram e{t) 



To test the validity of the universal power law behav- 
ior of as a function of the distance from the criti- 
cal load given by Eq. (11), simulations were performed 
with various disorder distributions, i.e. besides the uni- 
form distribution the Weibull distribution of the form 
P{e) = 1 — exp [— {e/\Y] was employed. The value of 
the characteristic strain A was set to one, and the shape 
of the distribution was controlled by varying the value of 
p. The results are presented in Fig. 3, where an excellent 
agreement of the simulations and the analytic results can 
be observed. Fig. 3 supports that the exponent oitf as 
a function of (Tq — ac is universal, it does not depend on 
the specific form of the disorder distribution. 

To study the finite size scaling of the time to failure 
tf a. uniform distribution was used for the failure thresh- 
olds. The value of the external load was fixed above 
(Tc and the number of fibers N was varied from 5 x 10^ 
to 10''. Averages were calculated over 10** samples for 
each system size N. The results obtained by simulations 
are presented in Fig. 4, where an excellent agreement of 
simulations and analytic results can be observed for four 
orders of magnitude in the system size N. 
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B. Slowly relcLxing fibers 

1. Analytic model 

Another important microscopic mechanism which can 
lead to macroscopic creep is the slow relaxation follow- 
ing fiber failnrc. In this case, the components of tlie solid 
are linearly elastic until they break, however, after break- 
ing they undergo a slow relaxation process, which can be 
caused, for instance, by the sliding of broken fibers with 
respect to the matrix material or by the creeeping ma- 
trix. To take into account this effect, our approach is 
based on the model introduced in Refs. [9,10], where the 
response of a viscoelastic-plastic matrix reinforced with 
elastic and also viscoelastic fibers have been studied. 



surrounding matrix material are modeled by Maxwell el- 
ements as illustrated in Fig. 5, i.e. they are conceived as 
a serial coupling of a spring and a dashpot which results 
in a non-linear response 



(16) 



where ct, and ey, denote the time dependent load and de- 
formation of a broken fiber, respectively. The relaxation 
of the broken fiber is characterized by three parameters 
Eb.B, and m, where Eb is the effective stiffness of a bro- 
ken fiber, and the exponent m characterizes the strength 
of nou-liuearity of the element. We study the behavior 
of the system for the region m > 1. 




FIG. 4. The size dependent lifetime of the bundle. The 
inset shows the lifetime of relatively small systems on a linear 
plot, while in the main figure the difference of the lifetime of 
finite bundles and the infinite one obtained from Eq. (14) can 
be seen on a double logarithmic plot. The slope of the fitted 
straight fine is -0.99 ± 0.02. 



The model consists of N parallel fibers, which break 
in a stress controlled way, i.e. subjecting a bundle to 
a constant external load fibers break during the time 
evolution of the system when the local load on them 
exceeds a stochastically distributed breaking threshold 
ai, i = 1, . . . ,N. Intact fibers are assumed to be lin- 
early elastic i.e. a = EfSf holds until they break, and 
hence, for the deformation rate it applies 



£f = 



Ef' 



(15) 



Here £/ denotes the strain and Ej is the Young modulus 
of intact fibers, respectively. The main assumption of the 
model is that when a fiber breaks its load does not drop 
to zero instantaneously, instead it undergoes a slow re- 
laxation process introducing a time scale into the system. 
In order to capture this effect, the broken fibers with the 





FIG. 5. The model solid when intact fibers are linearly elas- 
tic, and the broken ones with the surrrounding matrix are 
modeled by Maxwell elements. 



Assuming global load sharing for the load redistri- 
bution, the constitutive equation describing the macro- 
scopic elastic behavior of the composite reads as 



ao = ait) [1 - P(a(t))] + (Tb{t)P{a{t)). 



(17) 



Eq. (17) takes into account that broken fibers carry also 
a certain amount of load (7b{t), furthermore, P{a{t)) and 
1 — P{a{t)) denote the fraction of broken and intact fibers 
at time t, respectively [25]. It can be seen from Eq. (17) 
that under a constant external load ctq, the load of in- 
tact fibers a will also be time dependent due to the slow 
relaxation of the broken ones. 

Due to the boimdary condition illustrated in Fig. 5, 
the two time derivatives have to be always equal 



Ef = Eb- 



(18) 



The differential equation governing the time evolution of 

the system can be obtained by expressing ay, in terms 
of a from Eq. (17) and substituting it into Eq. (16) and 
finally into Eq. (18) 
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B 



1 1 

(To - (T [1 



1 



P{a) ^ P{ar 



(O- - (To) 



= (19) 



P{a) 



P{a) 



In order to determine the initial condition for the inte- 
gration of Eq. (19) the breaking process of fibers has to 
be analyzed. Subjecting the undamaged specimen to an 
external stress CTq all the fibers attain this stress value im- 
mediatelly due to the linear elastic response. Hence the 
time evolution of the system can be obtained by inte- 
grating Eq. (19) with the initial condition a{t = 0) = CTo- 
Since intact fibers are linearly elastic, the deformation- 
time history e{t) of the model can be deduced as e{t) = 
a{t)/Ef, which has an initial jump to £<, = ag/Ef. It 
follows that those fibers which have breaking thresholds 
(Tj smaller than the externally imposed CTq immediately 
break. 

To characterize the macroscopic behavior of the com- 
posite the solutions a{t) of Eq. (19) have to be analyzed 
at different values of the external load ctq. Similarly to 
the previous model, two different regimes of (7{t) can be 
distinguished depending on the value of CTq: if the ex- 
ternal load falls below a critical value CTc a stationary 
solution ag of the governing equation exists which can be 
obtained by setting ct = in Eq. (19) 



iTo = (Ts [1 - P((Js 



(20) 



This means that until Eq. (20) can be solved for as the 
solution (T{t) of Eq. (19) converges asymptotically to Us 
resulting in an infinite lifetime tf oi the composite. Note 
that Eq. (20) provides also the asymptotic constitutive 
behavior of the model which can be measured by qua- 
sistatic loading. If the external load falls above the criti- 
cal value the deformation rate i = &/Ef remains always 
positive resulting in a macroscopic rupture in a finite time 
tf. It follows from Eq. (20) that the critical load cTc of 
creep rupture coincides with the static fracture strength 
of the composite. 

The behavior of the system shows again universal as- 
pects in the vicinity of the critical point. Below the criti- 
cal point the relaxation of a{t) to the stationary solution 
as is governed by a differential equation of the form 



di 



(21) 



where 6 denotes the difference 5{t) = ag — <y{t). Hence, 
the characteristic time scale r of the relaxation process 
only emerges if m = 1, furthermore, in this case also 
T ~ (cTc — tTo)~^/^ holds when approaching the critical 
point. However, for m > 1 the relaxation process is char- 
acterized by d{t) = at^/^~™, where a — > with ao ac- 
Similarly to the previous model, it can also be shown 
that the lifetime i/ of the bundle has a power law di- 
vergence when the external load approaches the critical 
point from above 



tf ~ ((70 - (Tc) 



-(m-1/2) 



for (To > (Tc- 



(22) 



The exponent is universal in the sense that it is indepen- 
dent on the disorder distribution, however, it depends 
on the stress exponent m, which characterizes the non- 
linearity of broken fibers. 



2. Simulation technique 

Subjecting a finite bundle of TV fibers to ao external 
stress those fibers whose failure threshold falls below cTq 
break immediately. The number Ng of intially breaking 
fibers can be estimated from the disorder distribution as 
No ~ NP{ao). In the presence of broken fibers the sys- 
tem slows down and the remaining fibers of the bundle 
break one-by-one in the increasing order of their breaking 
thresholds a^a+i < <7No+2 < <7n- In order to construct 
an efficient simulation technique one has to determine the 
time elapsed between two consecutive breakings during 
the creep process. 

The macroscopic constitutive equation for a system of 
N fibers when i fibers have already failed can be written 
as 



N 



(23) 




tf 

_i \ 1 \ 

600 800 1000 1200 1400 
t/to 

FIG. 6. £ as a function of t for several value of Oo below 
and above Oc- N = lO' fibers were used. 



Making use of Eqs. (16,18), the differential equation 
describing the time evolution of the load of intact fibers 
(T can be cast in the form 



a 



J--^fl-^ 
Ef Eb\ i 



B 



(24) 



where fi{x) is introduced for brevity as 



6 



N 



(25) 



The time At((Ti, di+i) elapsed between the breaking of 
the ith and i + 1th fiber can be determined by integrat- 
ing Eq. (24) from cr, to CTj+i, which yields for m ^ 1 



K.. 



[m — 1) 

and the multiplication factor Ki reads as 



\1— ml 



N 



N - 



1 
B 



1 



1 



1 - 



N 



(26) 



(27) 



For TO = 1 the corresponding equation has the form 

At(<Ti, CTi+i) = Ki [ln/i(ai+i) - ln/i(<Ti)] . (28) 

Then the simulation proceeds as in the case of vis- 
coelastic bundles but in the above formulas the number 
of broken fibers i varies as z = No, A^o + 1, • • • , ^ — 1, so 
the time as a function of a can be obtained as t{(Ji) = 
'l2]=Na+i^^i^i'^i+'^) from which the deformation as a 
function of time Si{t) can be determined, since a = EfSf 
always holds. The lifetime tf of the system can be ob- 
tained by summing up all the At's. 




FIG. 7. Lifetime tj as a function of the distance from the 
critical point CTo — ffc for two different values of the parameter 
m. The number of fibers in the bundle was taken N = IQ' . 



For the purpose of explicit calculations a uniform dis- 
tribution was prescribed for the breaking thresholds cr, 
between and 1. The deformation as a function of time 
is plotted in Fig. 6 for several different values of the ex- 
ternal load below and above the critical load. Similarly 
to the previous model the two regimes of the creeping 
system can be clearly distinguished. 

To study the behavior of the time to failure as a func- 
tion of the distance from the critical point, simulations 



were performed for several different values of the expo- 
nent TO. In Fig. 7 the results are presented for to = 1.5 
and TO = 2.5. The slope of the fitted straight lines agrees 
very well with the analytic predictions of Eq. (22). 

The size scaling of the time to failure tf was ana- 
lyzed by simulating the creep rupture of bundles of size 
N = 5 ■ 10^ — 10^ setting a uniform distribution for the 
breaking thresholds. We found that tf{N) converges to 
the lifetime of the infinite system tf{oo) according to the 
universal law Eq. (14) independently on the value of the 
exponent m. In Fig. 8 the best fit was obtained for both 
curves with slope —1.0 ± 0.05 for both m values. 




N 

FIG. 8. Size dependence of the lifetime tf for two different 
values of the parameter m. 



III. CONCLUSIONS 

The creep rupture of fibrous materials occurring under 
a steady external load is microscopically a rather complex 
phenomenon depending on a diversity of possible mate- 
rial specific mechanisms. Therefore, on the one hand, it 
is impossible to work out a general theoretical framework 
which takes into account all the features of the process 
and has predictive power, and on the other hand, it is 
very important to reveal universal aspects of the creep 
process, which do not depend on specific material prop- 
erties relevant at the micro level. 

In the present paper we studied the creep rupture of 
fibrous materials in the framework of fiber bundle models 
taking into account two possible microscopic mechanisms 
of creep: (i) in the first approach the fibers themselves 
are viscoelastic and they break when their deformation 
exceeds a stochastically distributed threshold value, (m) 
in the second model the fibers are linearly elastic until 
they break, however, after breaking their relaxation is not 
instantaneous but the creeping matrix introduces an in- 
trinsic time scale for the relaxation. The first model can 
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be relevant for natural fiber composites like wood which 
are composed of viscoclastic fibers [f-4], while the sec- 
ond model can provide an adequate description of metal 
matrix composites reinforced by brittle fibers [29]. Ana- 
lytical and numerical calculations showed in both models 
that increasing the external load on a specimen a tran- 
sition takes place from a partially failed state of infinite 
lifetime to a state where global failure occurs at a finite 
time. The critical load turned to be the static fracture 
strength of the material. It was found that although 
the details of the two model systems are different, in the 
vicinity of the critical point they exhibit universal be- 
haviour which can be relevant for experiments. 
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V. APPENDIX 



where Xi is defined implicitly by P{xi) = i/{N + 1). Ap- 
plying Eq. (31) to Eq. (29) the resulting summation can 
be approximated by integrals replacing i /N by the equiv- 
alent P{xi){l + 1/N). Neglecting corrections higher order 
ml/N after straightforward calculations we arrive at 

{tf)^ j dxd2G{P{x),x) 

+ 2^ / dxP{x){l-P{x))dfd2G{P{x),x), (32) 

where d2G{y,x) = dG{y,x)/dx, and d1d2G{y,x) = 
d^G{y, x)/d^ydx. Substituting the actual form of G(y, x) 
for a specific model the complete form of the size scaling 
of lifetime can be obtained. However, it can be seen in 
the general expression Eq. (32) that the first term pro- 
vides the lifetime of the infinite bundle and the only size 
dependence is in the prefactor of the second term. Eq. 
(32) states that if the lifetime can be written in the form 
of Eq. (29) the lifetime of finite bundles converges to that 
of the infinite one asl/N with increasing number of fibers 
N. 



Here we provide the derivation of the average lifetime 
for the general case when the lifetime of a bundle with 
a specific realization of the disorder can be cast in the 
form 



N-l r 



G{j^,Xi+i)-G{j^,Xi) 



(29) 



i.e. tf is a sum of terms which depend on the number of 
broken fibers i and on a single breaking threshold .t, that 
can be given as strain or stress. The Xi-s are obtained 
by choosing N breaking thresholds independently from 
a cumulative probability distribution P{x) and putting 
them into increasing order. This treatment includes both 
models discussed in the present paper. The expectation 
value of a function f{xi) can be determined as 



I 



m 



{i-iy.{N-i)\ 

xp{x)f{x)dx. 



p{xy-' [i-p{x)] 



N- 



(30) 



The probability distribution in Eq. (30) that the value 
of the ith largest breaking threshold falls between x and 
X + dx has a sharp peak for large N values for each i. 
The above integration can be carried out by expanding 
the distribution about its peak. After expanding the re- 
sult in terms of 1/N and neglecting higher order terms 
we arrive at 

{f{Xi)) = f{Xi) + ^,P{Xi){l - P{Xi)) 



2N 
[P'ixiW 



f {Xi)P {Xi) 



[P'{x,W 



(31) 
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